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I. INTRODUCTION 



The idea of studying a gravitational system by replicating certain of its aspects in a laboratory environment through 
other, analogous, means has gained wide popularity today. One of the main motivations behind the current interest 
in these analog systems is to study the Hawking effect, which shows that a black hole evaporates through a process 
of spontaneous particle creation and emission [l| . It was later shown that this type of particle creation is not unique 
to black hole space-times Q but is a generic feature of any space-time possessing an apparent horizon Q. 

The process by which Hawking radiation is created is still not entirely understood, and seems virtually impossible 
to observe in gravitational systems. The search for a non-gravitational physical system that possesses some kind 
of analog to the Hawking effect originated with so-called "dumb holes" , comprising acoustic waves in moving fluids 
[3, HI, although a greater degree of success has been more recently achieved with surface waves, rather than acoustic 
waves, in moving fluids £6|— lS|] . See Refs. |9l. 1 1 01] for a review of other promising non-gravitational systems, including 
Bose-Einstein condensates and superfluids. It is alsopossible to consider electromagnetic analog systems such as 
"slow light" propagating in hi ghly dispersive media [T5L [Tlj , nonlinear electrodynamics [l2| , and light propagating 
in moving dielectric fluids (l3l Il4j. While a physical realization for electromagnetic analog space-times would pose 
challenges for reproducing the Hawking effect jl5lll6j|. it is - in principle - possible, and would certainly be useful for 
studying classical effects. 

A somewhat simpler electromagnetic approach is to identify the behavior of electromagnetic fields on a curved 
space-time manifold with the behavior of electromagnetic fields in a dielectric material residing in Minkowski space- 
time. This idea seems to have been first advanced by Eddington to argue that the deflection of light passing through 
the curved space-time near a massive object would be indistinguishable from the deflection of light passing by the 
same object if the space-time were Minkowskian but filled with an appropriate dielectric material [l7| . Gordon then 
examined the inverse problem of finding an equivalent curved space-time to describe a given dielectric material by 
an effective "optical" metric [l8|]. In a similar vein it was noticed by Plebanski that the form of the constitutive 
equations for electromagnetic fields in empty, curved, space-time is formally equivalent to the constitutive equations 
for electromagnetic fields in an appropriate dielectric This formal equivalence has been exploited to study 

gravitational systems in only a few examples [2(| HH • 

Two issues have contributed to the underdevelopment of this approach. One is that the materials required to mimic 
a gravitational system generally have somewhat unnatural characteristics, such as equal permeability and permittivity. 
Another is that, as cautioned by Plebanski himself 19], the Plebanski equations are not strictly covariant, which calls 
into question the limits of their applicability. However, the recent development of man made "metamaterials" [22l. |23| 
allows for wide ranging control of material parameters, to the extent of even achieving a negative refractive index 



24| . Here we address the second issue by describing a completely covariant approach to dielectric analog space-times 



that generalizes the Plebanski approach. Clearly there are close connections between metamaterials and analog space- 
times, and the method described here is closely related to a particula r appro ach [25l |26| in "transformation optics" , 
the emerging field of transformation-based metamaterial engineering }27H30l ]. While these are still the early days of 
metamaterial exploration, the approach has already shown the potential to address fundamental questions [3 II ] . 

This paper is organized as follows. Section[H]gives a brief review of electrodynamics in linear dielectric materials from 
a completely covariant, manifestly 4-dimensional standpoint. Section IIIII addresses the mapping of electromagnetic 
fields from a curved, vacuum, manifold to a material residing in Minkowski space-time. The characteristic parameters 
describing the analog material are determined as a function of position. Section HVl reconsiders some examples known 
from the literature, not only demonstrating agreement with results obtained from the Plebanski equations but also 
illustrating some additional mappings allowed by this new approach. In Sec. [V] we examine more carefully the 
correspondence between null geodesies in the vacuum space-time and trajectories of light in material. Section IVII 
considers the dielectric analog of the exterior Kerr geometry. We conclude with Sec. IVIII 



II. CLASSICAL ELECTRODYNAMICS 



The covariant description of electrodynamics in vacuum is thoroughly described in Ref. [32j |. whose notation and 
sign convention we follow, as well as a multitude of other, very readable, sources, such as Ref. (33j . Here we briefly 
summarize a covariant description of electrodynamics in linear dielectric media that has been recently given [26| , but 
for a more complete study of electrodynamics in dielectric media see Ref. [34j . 

Assume space-time to consist of a manifold M that possesses a metric g. The covariant approach to electrodynamics 
combines the electric field E and magnetic flux B into a single mathematical object, the field strength tensor F, that 
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in a locally Lorentz frame with Cartesian coordinates has components that can be represented as a matrix 
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Additionally, the electric flux D and magnetic field H are combined into the excitation tensor G, that can also be 
represented as a matrix 



_ -n x u u z -u y , , 






H x 




H z 


—H x 





D z 


-Dy 


-By 


-D z 





D x 


-H z 


Dy 


-D x 






The covariant Maxwell equations are expressed as 

dF = 0, dG = J, (3) 

where d is the exterior derivative, and J is the charge-current 3-form. 

Furthermore, in a linear dielectric medium there exists a relationship between F and G given by the constitutive 
equation (25l. l2q 

G = X(*F), (4) 

that in component form reads 

Gfiv = X^J*^ * a p P E a p- (5) 

In Eq. ([4}, * is the Hodge dual on (Af, g), which for present purposes is to be understood as a map from 2-forms to 
2-forms that has component form 



* 



"' -- /l 5|e^ ff V' / - (6) 
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The tensor x contains information on the dielectric material's properties (i.e. permittivity, permeability, and magneto- 
electric couplings), and can be thought of as representing an averaging over all the material contributions to an action 
that describes a more fundamental quantum field theory [35l l36j . We require % to be independently antisymmetric 
on its first two and last two indices, and in vacuum x(*F) = *F. This last condition means that the classical vacuum 
is treated as a linear dielectric with trivial \, recovering the usual constitutive relations in vacuum. 

The components of the constitutive equation provide a set of six independent equations that, in a local frame, can 
be collected in the form 

H a = (ji- l ) a b B b + (ix*) a b E b , D a = (e*) a b E b + (7 2 *) Q b S 6 . (7) 

where we use the notation a to denote a 3 x 3 matrix. Rearranging these to 

B a = (fi) a b B b + (7i) a b ^ b , D a = (s) a »E b + (i 2 ) a b B b . (8) 

gives a representation that may be more familiar. These three-dimensional representations of the completely covariant 
Eq. ([¥]) are essentially equivalent, and it is a simple matter to switch between them using the relations 

e = e* - j2*Mi*> 7i = ~Mi*> 7~2 = 72*A- (9) 



However, one should be aware that these 3x3 matrices are not tensors but simply components of x that have been 
collected into matrices. By matching with the usual representations above, one finds the matrix representation of x 
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where the * indicates entries that are antisymmetric on either the first or second set of indices on x 1 s TP ■ Equation 
(fTOf represents x as a matrix of matrices; the first two indices of x 1 ^ p give the component of the large matrix, 
which is itself a matrix described by the second set of indices on x- 



III. DIELECTRIC ANALOG SPACE-TIMES 



The idea behind a dielectric analog space-time is to identify the behavior of electromagnetic fields on a curved space- 
time manifold with the behavior of electromagnetic fields in a dielectric material residing in Minkowski space-time. 
More precisely, let (M, g) be a curved, vacuum, space-time manifold and let F and G be electromagnetic fields in the 
space-time. In particular, it is assumed that g is a solution of Einstein's equations, rather than the Einstein-Maxwell 
equations, so g describes the background space-time onto which some electromagnetic fields have been placed. Thus 
we want to identify the null geodesies of the curved, vacuum, space-time with trajectories through the equivalent 
dielectric medium as determined by the wave equation in that medium. Because the curved space-time is vacuum, 

X = \vac- 

The question that we would like to answer is: given a curved space-time, how do we determine the dielectric 
material that mimics this curved space-time? The answer must come in the form of some characteristic material 
parameters such as permeability and permittivity. To see how to obtain such material parameters as a function of 
position, consider a map T : M —> M relating the electromagnetic fields in the material with electromagnetic fields 
in the curved space-time, as depicted in Fig. [TJ 



(M,g,*) 

;f,g, x ) g = x (*f; 

Minkowski with materials 



A A 



T 



;f,g, x ) g = x(*f; 

Curved vacuum manifold 



FIG. 1. Electromagnetic fields F and G are pulled back from a patch of a curved, vacuum manifold to a material residing in 
Minkowski space-time by using the pullback T* . 

Because the metrics g and g are known in each space-time, we demand that T acts only on the electromagnetic 
fields via its pullback, T*, by 

T* (if) = F. (11) 
Consider now the pulled-back fields. At a point x G M, 

G x = T* (G Tix) ) = T* (xn,) ° *T(x) ° Fr W ) ■ (12) 
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But since G = x(*F) at x £ M, it is also true that 



(13) 



It follows that the right hand sides of Eqs. ([12"]) and (fT3")) must be equal. To clearly see what is going on, consider the 
action of G x on a bi- vector V x £ T%(M), then 



Xx o + x o t* If 
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(14) 



Let be the Jacobian matrix of T, which is the matrix representation of dT. Then Eq. (fl"4|) can be written in 
component form as [25j 
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where we have explicitly indicated where each object is evaluated. Eliminating F and ~V X from both sides and solving 
for x as a function of a: S M gives 



a/3 



T{x) 



f.LV 



T(x) 



p *n6 



(16) 



Here A -1 denotes the matrix inverse of A, both A and A -1 are evaluated at x, and we have used the fact that on a 
4-dimensional Lorentzian manifold, acting twice with * returns the negative, **F = — F. If the initial space-time is 
vacuum, then = & and 
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Equation (I17p gives the permeability, permittivity, and magneto-electric couplings of a material that mimics a 
curved, vacuum, manifold. The metric g of the curved space-time is arbitrary; in particular g need not be static Q. 
If the map 7" is chosen to be the trivial coordinate identification 7o (x) = x' (technically this map is not the identity, 
since the manifolds are different), then 



(18) 



This choice of T recovers results obtained with the Plebanski equations, but in general T need not be the trivial map 
7o, allowing for greater freedom in identifying materials. For example, we show below how the exterior Schwarzschild 
geometry may be mapped into a spherical shell of finite thickness, albeit at the expense of introducing divergences in 
the material parameters. It should be pointed out that the space-time represented by some dielectric configuration 
is not unique. This is because the Hodge dual ★ is invariant under conformal transformation of the metric g — > fi 2 g. 
Therefore a dielectric analog space-time actually represents an equivalence class of conformally related space-times. 



IV. RINDLER WEDGE AND EXTERIOR SCHWARZSCHILD GEOMETRY 



A. Rindler Wedge 

As a simple example, consider the Rindler geometry described by the metric 

ds 2 = -a 2 (x') 2 (dt') 2 + (dx') 2 + {dy'f + (dz 1 ) 2 . (19) 

The regions x' < and x' > correspond to the left and right Rindler wedges, respectively. In this and what follows, 
primed coordinates denote the curved, vacuum manifold, while unprimed coordinates denote those of Minkowski 
space-time with dielectric media. Thus the desired material parameters will be described as functions of unprimed 
coordinates. Choosing T = 7o, Eq. ([18} returns a material described by scalar- valued permeability and permittivity 

£ = u = — (20) 
ax 

in agreement with results obtained from the Plebanski equations [2l[ . Notice that the material parameters diverge at 
the singular surface x = 0. 
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B. Exterior Schwarzschild Geometry 

The dielectric analog of spherically symmetric gravitational systems in isotropic coordinates was first studied by 
De Felice [20], while Reznik 21] mentions the material parameters representing the exterior Schwarzschild geometry. 
Taking the line element in Schwarzschild coordinates 

ds 2 =-(l-^\ {dt 'f + (i _ *^f\ 1 {dr >f + (r'f(d0') 2 + (r'f sin 2 6'(d V 'f (21) 

and the map %(t, r, 9, ip) — (f , r', 9', <p') — (t, r, 9, ip) returns the material parameters 

/ 2A/\ _1 

E rr = fl rr = 1, Egg = E vv — flgg — fl vv = I 1 — J (22) 

in the representation of Eq. (|5|). Once again the horizon is characterized by divergent material parameters. Also note 
that in both of the previous examples the material fills the entire space. However, by adopting some map T ^ 7o, we 
may gain some freedom in specifying the boundaries of the material. For example, adopting 

T(t, r, 9, V ) = {if, r',9', ip') = (t, 2A fj r a) , 9, <pj (23) 

maps the exterior Schwarzschild geometry onto a spherical shell a < r < b. Using this map in Eq. (|17|) now leads to 

_2M(b~a) _ 2M(b-a) 2 

E rr — [l rr — 5 , Egg — £ w — flgg — — 7 rr; (24) 

r z (r — a)(b — r) z 

In a sense this result is worse in terms of the material, as its permeability and permittivity now diverge at both the 
inner and outer surfaces of the shell. However, this behavior is not unexpected, and is similar to the example of 
compactifying a flat space-time containing an electric charge, which results in the appearance of an image charge on 
the boundary. 

We may go a step further and try to remove the divergences in e — fi for the compactified exterior Schwarzschild 
analog by modulating the time coordinate with a map such as 

n.r,>M - e-y,«v> - ( ,"'■'-•'' gM,,, y ) . ,25, 

\(r — a)(b — r) z b — r J 

We expect that this kind of map should eliminate the infinite blue shift by redefining time as the wavelength decreases 
towards the singular surfaces. Turning the crank on Eq. (|17p returns 

(r — a)(b — r) 2 

E rr = fX rr = — — , Egg — £m<p = = l^vv = 1; (26) 

(o — a)r z 

which is finite, as desired. However, this type of spatially dependent time transformation generically introduces 
magneto-electric couplings [251 ] . which in this case turn out to be 



Thus the divergences have been shifted from the permeability and permittivity to divergences in new magneto-electric 
couplings. In the black hole space-time, ingoing light rays will be lost behind the horizon, so it is not surprising that 
trying to mimic this behavior in a material should require these essential singularities to appear somewhere. Note that 
not only do the magneto-electric couplings diverge at the inner and outer surface of the material shell, but there is 
also a sign change when 71, and 72 momentarily vanish at r = 2a + b . On the Schwarzschild manifold this corresponds 
to the point r' = 3M, which is precisely the radius of the photon sphere in the Schwarzschild geometry. 

V. WAVE EQUATION AND GEOMETRIC OPTICS IN MATERIALS 

Before considering any further examples, return first to the constitutive equations and material tensor X- As they 
stand, Eqs. ([7]) and © are somewhat misleading, because while E a and H a are components of a 1-form, D a and B a 
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are really selected components of the 2-forms D ab dx a A dx b and B ab dx a A dx b . Thus, a more appropriate version of 
Eq. (0 would be something like 



(*vB) a = (fi) a b H b + (i 1 ) a b E b 



(28a) 



(* E D) Q = {e) b E b + (i 2 ) a b H b 



(28b) 



where is the Hodge dual on an appropriate 3-dimensional spatial hypersurface. But this requires that we resolve the 
space-time into space and time components, selecting an observer to define a direction of time. The spatial hypersurface 
is then orthogonal to the selected direction of time. Transforming to the local frame of the selected observer we can 
make the identifications of Eq. (|10[) and then give the constitutive equations a 3-dimensional representation. 

In the Schwarzschild example above we expressed the results in spherical coordinates in the representation of Eq. 
(|8"|), but in dealing with the wave equation and determining the trajectories of light through a material it is the 
covariant x that is of real importance. To see this, return again to Maxwell's inhomogeneous equation in the absence 
of sources 



dG = 0. 



(29) 



Writing the exterior derivative of a 2-form in terms of the codiffcrcntial 6 — — * d-k, and using both the constitutive 
relation Eq. Q and the definition F = dA, where A is the 4-potential, it follows that 



S (*x * dA) = 

is the covariant wave equation for the 4-potential [37] . We may reassure ourselves that since Xv 
we recover the usual vacuum wave equation in the Lorenz gauge 8 A. = 0, 

MA = AA = 0, 



(30) 



★ m vacuum, 



(31) 



where A is the Laplace-De Rham operator. To see that the wave equation in the vacuum Schwarzschild space-time is 
equivalent to the wave equation in the analog material residing in Minkowski space-time, let us be more explicit with 
the notation and denote *s and *m as the Hodge dual in the Schwarzschild and Minkowski space-times respectively. 
For illustrative purposes let us further assume the map 7o- Then from Eqs. ([29)) and (|T8|) it follows that 



d(**A/F) = d(x V ac *sF), 



(32) 



demonstrating the formal equivalence of the analog model on the level of the wave equation. Thus, it is sufficient to 
describe the components of which may then be transformed to the local frame of any observer. 

We may go a step further in examining the covariant wave equation, and demonstrate that in the limit of geometric 
optics the trajectories of light r ays in the equivalent analog media are equivalent to the null geodesies of the vacuum 
space-time (see also e.g., Refs. [2(J H3])- For this purpose let us assume the geometric optics ansatz for A, 



— A^t 



(33) 



where e is some dimensionless parameter. Using this ansatz, Eq. (|29[) may be calculated explicitly, and the geometric 
optics limit is obtained by taking e — > [Hf . We do this for the analog exterior Schwarzschild geometry described by 
Eq. (j22j) . Consider first those modes with fc M = (kt,k r , 0, 0), for which Eq. ([29]) returns 
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Thus we find the condition on fc M that 



, 2M\ , 
kt = ( 1 — ) Kr, 



(35) 



which is precisely the condition determined by the requirement that k^k^ = for radial null geodesies in the vacuum 
Schwarzschild space-time. For a more general k^ — (kt, k r , kg, k v ), calculating Eq. (|29[) results in a more complicated 
expression that can be written as a system of four equations 
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(36a) 



8 



Ath [l-'m + A r k r (1 — 



A h 

r 2 sin 2 i 



A, 



;i-^r 1 + fc. 2 (i-2M) + 



hp 



r 2 sin 2 i 



= 0, (36b) 



-A t k t (1 



2M1 



Aeke ^ A v k v 
r 2 r 2 sin 2 £ 



-k 2 t (i 



2M^ 



r 2 sin 9 



0, 



(36c) 



2i\/\ 



+ 



1 - 



2M 



4 £• 
r 2 sin 2 6 



the simultaneous solution of which requires 
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(36d) 



(37) 



corresponding to the condition k^k^ — for a general null geodesic in the Schwarzschild space-time. We have thus 
demonstrated that in the limit of geometric optics, the dielectric analog Schwarzschild black hole accurately mimics 
the behavior of light propagation in the vacuum space-time. 



VI. EXTERIOR KERR GEOMETRY 



Consider next the exterior Kerr geometry in Boyer-Lindquist coordinates described by the line element 

2Mr'\,,, N9 AMar 1 .,„,,,,, p 2 , , „ „,„ M S 



where 
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(38) 



(39) 



and a is the ratio of angular momentum to mass, so that L — aM . A naive calculation using the Plebanski equations 
[l9l . l20l | , or using Eq. (jl8]l and the representation Eq. (jHJ , returns material parameters like 



2Mr' 



(40) 



which diverge on the time-like limit surface (TLS), the boundary of the ergoregion. This is a curious result. The Boyer- 
Lindquist coordinates are well behaved at the TLS, so we should not expect divergent behavior at the corresponding 
material points. However, in the previous section we have seen that x is really the only meaningful covariant quantity, 
and the map 7o provides the analog Kerr parameter in spherical coordinates as 
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(41) 



where a * denotes entries that are antisymmetric on either the first or second pair of indices of X 1 $' P ■ It is clear that 
X is completely regular and finite across the TLS, but that some of the components of x change sign at the TLS. 
When the inverse of these components are used in the representation of Eq. ((5]) a divergence appears. In terms of the 
propagation of light, which depends on x, everything is well behaved at the TLS. The only divergences occur on the 
horizon, A = 0, as expected. 
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VII. CONCLUSIONS 



Using a manifestly 4-dimensional covariant description of electrodynamics, a natural and covariant method of iden- 
tifying a dielectric analog of a curved space-time has been developed. First, the quantities that describe macroscopic 
electrodynamics in materials were encoded in such a way that the vacuum may be considered as just another dielectric 
material. Then, by pulling back the fields from the curved space-time to Minkowski space-time with a material, the 
requisite material parameters could be solved for. However, it was found that the map T, whose pullback T* operates 
on the fields, is not uniquely defined. 

The non-uniqueness of T is somewhat unsettling, and it would be nice if there were a natural choice for T. 
On the other hand, the greater freedom in selecting T allows for greater manipulation of the resultant material, 
where for example, it is possible to map the exterior Schwarzschild geometry to a finite spherical shell residing in 
Minkowski space-time. Such non-trivial mappings come at a price, such as more complicated material parameters, 
extra divergences, and induced magneto-electric couplings. 

In analyzing the analog exterior Kerr geometry, it was found that the "usual" representation of the constitutive 
relations requires careful consideration as it may lead to spurious divergences. However, the quantity with real 
relevance for the behavior of light in the material is x an( i the spurious divergences of the usual, 3-dimensional 
representation occur because they rely on the inverse of elements of x that may vanish at a point. 
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